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Abstract 

We report on the 2"5 (400 pc) resolution CO (J 
= 1 —> 0) observations covering the whole length of 
the bar in the strongly barred late-type spiral galaxy 
NGC 7479. CO emission is detected only along a 
dust lane that traverses the whole length of the bar, 
including the nucleus. The emission is strongest in 
the nucleus. The distribution of emission is clumpy 
along the bar outside the nucleus, and consists of 
gas complexes that are unlikely to be gravitationally 
bound. The CO kinematics within the bar consist 
of two separate components. A kinematically dis¬ 
tinct circumnuclear disk, < 500 pc in diameter, is un¬ 
dergoing predominantly circular motion with a max¬ 
imum rotational velocity of 245 km s“^ at a radius 
of 1" (160 pc). The CO-emitting gas in the bar out¬ 
side the circumnuclear disk has substantial noncir¬ 
cular motions which are consistent with a large ra¬ 
dial velocity component, directed inwards. The CO 
emission has a large velocity gradient across the bar 
dust lane, ranging from 0.5 to 1.9 km s“^ pc“^ after 
correcting for inclination and the projected velocity 
change across the dust lane is as high as 200 km s“^. 
This sharp velocity gradient is consistent with a shock 
front at the location of the bar dust lane. A compar¬ 
ison of Ha and CO kinematics across the dust lane 
shows that although the Ha emission is often observed 
both upstream and downstream from the dust lane, 
the CO emission is observed only where the velocity 
gradient is large. We also compare the observations 
with hydrodynamic models and discuss star formation 
along the bar. 

Subject headings: galaxies: evolution — galaxies: individ¬ 
ual (NGC 7479) — galaxies: ISM — galaxies: kinematics 


and dynamics — galaxies: starburst — galaxies: struc¬ 
ture 


1. Introduction 


To investigate the various factors that determine 
the molecular gas morphology and kinematics in the 
bars of spiral galaxies we need to obtain observations 
of galaxies which have molecular gas in both the nu¬ 
cleus and the bar. NGC 7479 is one of the few galaxies 
which have abundant molecular gas along most of the 
bar and in the nucleus ( [Quillen et al. 199^ ). It is an 
isolated, large, strongly barred SBc galaxy at a mod¬ 
erate distance (32 Mpc, assuming Hq = 75 km s“^, a 
heliocentric velocity ofV = 2371 km s“^ and no cor¬ 
rection for the local group velocities; 1" = 160 pc). It 
has both abundant atomic gas ( Laine fc Gottesman] 
I998|) and molecular gas (Young & Devereux 1991), 


and it has been classified as a starburst galaxy based 
on its large and centrally concentrated 10 /im flux 
(Devereux 1989). 

We have obtained new, high resolution (synthe¬ 
sized beam = 2'.'7x2'.'l, position angle = —85°) 2.6 
mm CO (J = 1 —> 0) observations of NGC 7479, 
and investigate the gas kinematics and distribution 
in the unusually gas-rich bar. We have detected a 
strong central molecular gas component and weaker 
emission along the bar, closely following the leading 
dust lane, as previo usly shown by the lo wer resolu¬ 
tion (7") CO maps ( [Quillen et al. 19% ). Our new 
data have greater sensitivity and spatial resolution. 
The emphasis in this paper is on the interesting CO 
kinematics in both the nuclear area and in the off¬ 
set gas/dust lane along the bar. Other papers in this 
series address the distribution and kinematics of the 
neutral, atomic hydrogen gas in NGC 7479 (Laine 
fc Gottesman 1998 ), the stellar bar pattern speed of 
NGC 7479 (Laine, Shlosman, & Heller 1998), and the 
minor merger model for NGC 7479 ([Laine & Hellei 

i^) . 


2. Observations 

The new CO observations were made at the Owens 
Valley Radio Observatory (OVRO) Millimeter Ar- 
ray.[| NGC 7479 was observed in three different con¬ 
figurations of the array. The observations were made 
on 1994 October 7, November 23-24, and 1995 Jan¬ 
uary 18-19. The total integration times at each of the 
three observed positions and other observing param¬ 
eters are given in Table ^ The total single sideband 
temperatures, including the effects of the atmosphere, 
were typically 350-500 K. Variations in the sky opac¬ 
ity and receiver gain were corrected through measure- 

^The Owens Valley Radio Observatory Millimeter array is op¬ 
erated by the California Institute of Technology with support 
from the National Science Foundation. 
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ments of an ambient temperature chopper wheel. The 
pointing was checked at the beginning of each track, 
and the strong quasar 3C 454.3 was observed every 19 
minutes to serve as a bandpass and gain calibrator, 
as well as a secondary flux calibrator. 

The calibration followed standard procedures. Af¬ 


ter I correcting for the atmospheric and instrumental 

gain variations, we produced channel and integrated 
maps of the CO emission. For mapping the visibility 
data, natural weighting was used to achieve the high¬ 
est sensitivity. We used 15 mJy beam“^ (1.5cr) as 
the stopping criterion for the CLEAN iterations. The 
map parameters are given in Table ^ We combined 
the three observed positions into a mosaic image, cor¬ 
recting for the attenuation of the primary beam. The 
velocity-smoothed channel maps of the CO emission 
are shown in Figure The measured line of sight 
CO velocities quoted in this paper are with respect 
to the LSR unless otherwise stated. The mapped re¬ 
gions have been drawn with great circles in Figure |^, 
which shows the integrated flux map. 

We also present data from Ha Fabry-Perot ob¬ 
servations, kindly given to us by Drs. Stuart Vo¬ 
gel and Michael Regan. These observations were 
made with the Maryland-Caltech Fabry-Perot Spec¬ 
trometer (see, e.g. Vogel et al. 1995) attached to 
the Cassegrain focus of the 1.5 m telescope at Palo- 
mar Observatory. The data were obtained on 1994 
September 29-30. 40 exposures were taken, each with 
a 500 second integration time and a pixel scale of 
1'.'88 pixel“^. To improve the signal to noise ratio, 
these data have been smoothed to 3'.'6 spatial res¬ 
olution. The velocity planes are separated by 12.1 
km s“^. The velocity uncertainty (rms) is estimated 
as 1-2 km s“^. The data reduction procedure fol¬ 
lowed that of Regan et al. (1996), who give a more 
detailed description. 

3. Results 

3.1. CO Flux and Distribution 

The integrated ^^CO (1-0) flux map of the whole 
mapped area and the mean velocity field of the cen¬ 
tral gas complex are displayed in Figures || and 
respectively. The unit of the integrated flux map has 
been converted to V(H 2 ) cm~^ with th e “standard 
conversion factor” (Bloemen et al. 1986 ; ^coville et 
d. 1987; Kenney fc Young 1989), 


V(H 2 ) cm-2 = 3.0 X 10^° J TB(CO)dv K km s"^ 

( 1 ) 

The CO flux spectrum is plotted in Figure The 
total integrated flux is 266±40 Jy km s“^, where the 
stated uncertainty is the rms value in the measured 
fluxes. In addition, systematic uncertainties in the 


flux of the secondary calibrator (15%) can lead to a 
systematic over- or underestimation of the fluxes. For 
comparison, the total, integrated single dish CO flux 
for the whole galaxy (4'x3') is 1050±200 Jy km s“^ 
( Young et al. 199^ ), but for the central 45" only, we 
estimate a CO flux of 430±65 Jy km s“^ (Kenney 

T^ ). 


At least 62% ± 20% of the flux in the central zone 
of the galaxy has been detected by our high resolution 
OVRO interferometer observations. Our interferom¬ 
eter observations have some sensitivity to structures 
up to 50" in size, although the sensitivity starts to 
decrease for structures with sizes larger than half the 
field of view (structures with sizes > 30"). The “miss¬ 
ing” flux could be associated with large scale struc¬ 
tures (> 30") for which our interferometer measure¬ 
ments are not fully sensitive. However, it is unlikely 
that a single velocity channel (width = 10 km s“^) 
would have continuous CO emission structures this 
large. More likely, our brightness temperature sensi¬ 
tivity (rms = 0.2 K, corresponding to an H 2 column 
density of about 10^^ molecules cm“^ averaged over 
the beam, using the standard conversion factor) is 
not good enough to detect the emission from regions 
with low beam-averaged CO surface densities. We 
will keep these limitations in mind when discussing 
the implications of our observations in § 4. 

Fig. I reveals that all the detected CO emission 
comes from a nearly linear ridge which is offset to¬ 
wards the leading edges of the bar (assuming the spi¬ 
ral arms are trailing). This ridge has a full-width- 
half-maximum (FWHM) value similar to the width of 
the beam (2'.'7), indicating that the gas/dust lane is 
unresolved. The dust lane, as measured from Ol'S res¬ 
olution optical and near-infrared images, has a width 
of about 1'.'5. Therefore, the high density molecu¬ 
lar gas ridge does not have a width larger than that 
of the dust lane. The CO emission intensity has 
a strong peak in the nucleus, but overall the emis¬ 
sion crosses the nucleus smoothly along position angle 
(PA) —10°, connecting the two linear, offset parts of 
the bar gas/dust ridge. 

Applying the standard CO-to-H 2 conversion fac¬ 
tor and multiplying the resulting masses by 1.35 to 
account for helium, the molecular mass within the 
central 1" (160 pc) radius is 3.3 x 10® Mq and 
the mean gas mass surface density in this area is 
4.1 X 10® Mq pc“^. This should be compared with 
the galactic center molecular gas mass of 3 x 10^ Mq 
within the central 500 pc (Dahmen et al. 1998) 
that has an average molecular gas surface density 
of 150 Mq pc“^. The corresponding hydrogen col¬ 
umn density in the nucleus of NGC 7479 is about 
1.6 X 10^® cm“^ which implies a beam-smoothed op¬ 
tical extinction (Ay) value of almost a hundred mag¬ 
nitudes. Extinction will be further discussed when 
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Table 1: Observing parameters at OVRO. 


Parameter 

Value 

Pointing center 1 

R.A. 23*' 02“ 25!99 

Dec. 12° 02' 24'.'2 

Pointing center 2 

R.A. 23^ 02“ 26!30 

Dec. 12° 03' 09'.'0 

Pointing center 3 

R.A. 23^ 02“ 26!61 

Dec. 12° 03' 53'.'7 

Central frequency 

114.37 GHz 

Central velocity (LSR) 

2340 km s“^ 

Configurations used 

A, B, C 

Total number of antennas 

A=5, B=6, C=6 

Field of view 

63" 

Spatial scales sampled^ 

2'.'5 - 50" 

Baseline lengths 

4.0 - 92.7 kA 

Integration time (min), p.c. 1 

350 

p.c. 2 

360 

p.c. 3 

360 

Absolute flux calibrator 

Neptune, Tb = 114 K 

Secondary flux calibrator 

3C454.3 

Velocity resolution 

10.4 km s“^ 

Total observed bandwidth 

464 MHz = 1206.4 km s"^ 


""See text in § 3.1 for more information. 


Table 2: Parameters of the channel maps. 


Parameter 

Value 

Size of synthesized beam 

2" 7 X 2'.'1 

Size of beam in parsecs®' 

262 X 204 

PA of synthesized beam 

00 

O 

Weighting used to make maps 

Natural 

Rms noise in channel maps 

0.011 Jy beam-i=0.19 K 

Theor. rms noise in channel maps 

0.0115 Jy beam-i=0.20 K 

Peak brightness temperature 

3.83 K 

Radio continuum flux density'^ 

< 4 mjy 


‘^Assuming a distance of 32 Mpc. 

5cr upper limit is given for a point source. 
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Fig. 1.— The 20 km s“^ channel maps of the naturally weighted CO cube. The contour levels are at (—3, 3, 5, 8, 
12, 15) X 0.021 Jy/beam. The star marks the AT-band nucleus. The velocity at the center of each channel is given 
in the upper left of each frame and the beam (2'.'7 x 2'.'1) at the bottom right of the last frame. Note the strong 
emission near the nucleus in most of the channels, due to the circularly rotating circumnuclear disk. 
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Fig. 2.— The naturally weighted integrated map of 
CO (1-0) emission. The contour and gray scale lev¬ 
els are (4, 8, 16, 32, 48, 96, 144) x 10^^ H 2 molecules 
cm“^. The contour around the CO emission is the 
17.1 mag arcsec“^ level of a iF-band image (Quillen 
et al. 1995), showing the location of the stellar bar. 
The PA of the stellar bar and that of the gaseous bar 
are different, since the gas is in a leading dust lane. 
The beam is shown at the bottom right corner of the 
map. The main CO clumps have been numbered. The 
beams around the three observed positions have been 
drawn with large circles. 



as** 2™26.5'' 26.O'" 

RA 


Fig. 3.— A magnified image of the central CO ve¬ 
locity field superposed on the integrated CO emission 
map. The CO emission peak is marked with a star in 
the nucleus. The gray scale levels have been chosen 
to show only the brightest CO emission. The isove¬ 
locity contours in this map do not accurately reflect 
the interesting gas kinematics in some regions. The 
Z-shaped contours near the center occur where the 
inwardly flowing bar gas meets the kinematically dis¬ 
tinct circumnuclear disk. These two distinct compo¬ 
nents are better displayed in Fig. Q. The large veloc¬ 
ity gradient across the bar dust lane, which we have 
marginally resolved, is better displayed in Fig. ||. 
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Fig. 4.— The spectrum of the CO emission detected 
with the OVRO millimeter interferometer. The chan¬ 
nel fluxes have uncertainties of the order of 5%. 
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addressing the star formation rates in § 4.2. The 
molecular gas mass in the central 7" diameter region 
is about 2 x 10® Mq and the total molecular mass de¬ 
rived from our observations is 4 x 10® Mq. Roughly 
half (2 X 10® Mq) of the detected flux (and mass) lies 
along the bar outside the central 7" diameter region. 

The detected CO emission along the bar is highly 
clumpy. In the integrated flux map (Fig. ||) at least 
six separate clumps in addition to the nuclear gas 
complex can be identified. The masses of the clumps 
along the bar range from 9 x 10® Mq to 2x 10® Mq. 
However, the use of the standard CO-to-H 2 conver¬ 
sion factor may not be justified in strongly shocked 
gas along the bar. The temperature and density are 
likely to be higher in the shock, although these effects 
tend to offset each other. If the clouds in the shock 
are not gravitationally bound (see § 4.1.2) their large 
line widths cause the standard CO-to-H 2 conversion 
to overestimate the H 2 mass. Shock chemistry can 
also change the conversion factor. The lower limit 
for the conversion factor is obtained for optically thin 
CO emission which has a conversion factor 20 times 
smaller than the standard factor (cf. Bryant & Scov- 
ille 1996). Overall, it is difficult to ascertain how the 
conversion factor would change in bar shocks. 


The sizes of the clumps along the bar are too 
large to be individual self-gravitating giant molecu¬ 
lar clouds but their sizes and masses are similar to 
those of giant molecular associations (GMAs) found 
in the nearby spiral galaxies M 51, M 100 (Rand 
1993a, 1993b, 1995), and NGC 4414 ( [Sakamoto 19“^ . 
However, the clumps in NGC 7479 are kinematically 
very different from those observed in the other galax¬ 
ies (§ 4.1.2). 

The azimuthally averaged radial profile of the molec¬ 
ular gas is compared to the radial profiles of R- and 
AT-band light, Ha, H I, and A = 21 cm radio con¬ 
tinuum emission in Figure |[ The optical and near- 
infrared images were first smoothed to the spatial res¬ 
olution of our CO data (2'.'5). Ho wever, the radio con- 
tinuum and the H I data (from [Laine fc Gottesman 
19S(^ have resolutions of 4" and 15", respectively. 


Within the central 7" the CO emission intensity 
rises by more than an order of magnitude, as in 
NGC 3504 (Kenney, Carlstrom, & Young 1993). Un¬ 
like NGC 3504, the R- and AT-band profiles are shal¬ 
lower than the CO profile in this area. This may 
reflect the lack of an intense starburst in the nucleus 
of NGC 7479, its relatively small bulge component 
or large extinction. The CO to Ha intensity ratio 
is approximately constant in the central 2", where 
the nucl ear gas disk lies, as in NGC 3504 ( Kenney et 
al. 1992). Outside this inner region NGC 3504 and 
NGC 7479 show opposite behaviours: in NGC 3504 
the Ha gradient becomes steeper than the CO inten¬ 
sity profile, whereas in NGC 7479 the opposite is true. 


This could mean that the star formation efficiency in¬ 
creases outwards from the center in NGC 7479. Un¬ 
fortunately, there are no suitable radio continuum 
measurements which would allow us to assess the ef¬ 
fects of extinction. In other galaxies comparisons be¬ 
tween the 6 cm thermal emission and Ha have shown 
that extinction can change the “picture” (e.g. M 100; 
Garcfa-Burillo et al. 1998). 

The radio continuum emission is very strongly 
peaked in the nucleus. These data are shown at 
4" resolution, so the gradient at a higher resolution 
is likely to be even larger. If the radio continuum is 
regarded as a star formation tracer, the star forma¬ 
tion efficiency decreases with radius as in NGC 3504. 
However, the nucleus of NGC 7479 is slightly active 
(LINER), and the radio continuum may also have a 
contribution from the nucleus, not related to star for¬ 
mation. 

Figure ||6 shows the striking difference in the dis¬ 
tribution of H I and CO emissions. It should be noted 
that our CO observations did not extend very far out 
into the spiral arms of NGC 7479. The atomic gas ap¬ 
pears to have been converted into the molecular form 
at the higher pressures of the central potential well, 
aided by the self-gravity of the gas. 


3.2. Gas motions 


The line of sight velocity field (Fig. ^ is complex 
and difficult to interpret. The contours in the nuclear 
few arcsec region are aligned perpendicular to the PA 
of the outer disk (22°; Laine & Gottesman 1998), con¬ 
sistent with solid body rotation. Beyond 3" radius 
in the south, the line of sight velocity contours are 
aligned roughly along the bar. This morphology of 
the velocity field is consistent with a large radial ve¬ 


locity component (along the bar; e.g., Lindblad et al 
1996|). The northern part of the bar does not have 


a sufficiently extended emission region in which we 
could estimate the character of the gas motion. 

The central 4" region of the velocity field where 
the circularly rotating gas and the gas flowing along 
the bar meet has a Z-shape. To better understand the 
gas motions in the bar we have made position-velocity 
plots at several different PAs and spatial locations. 

Position-velocity plots close to the major axis (at 
PA 25°) and along the bar are presented in Figures || 
and ^ respectively. Both plots are centered at the 
nuclear GO emission peak (R.A. 23® 02™ 26237, Dec. 
12° 03' ll'.'l; B1950.0) as determined from the map of 
integrated GO emission (Fig. ||), and represent planes 
(width 0'.'5) of the original data cube. Outside the 
central 2" region, the mean GO velocity within each 
emission region along the bar does not deviate from 
the mean velocity (2360 km s“^) of the nuclear region 
by more than 50 km s“^ (Fig. I). Fig. 0 also has a line 
which indicates the location where emission would 
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Fig. 5.— (a). The azimuthally averaged radial profiles of the CO emission, Ha emission, i?-band and iF-band light, 
and A = 21 cm radio continuum (RC) in the central area of NGC 7479. The data have been smoothed to the same 
spatial resolution of 275, except the radio continuum (4"). For the CO, the vertical axis gives the total molecular 
mass surface density in Mq pc“^, calculated using the standard CO-H 2 conversion factor. The other profiles, which 
have arbitrary units, have been normalized to agree at 5.0". The Ha data are from Young & Devereux (1991), the 
i?-band data were kindly provided by P. Martin and D. Friedli, taken at the Steward Observatory 2.3 m telescope 
on Kitt Peak. The RT-band data are from Quillen et al. (1995). The 21 cm radio continuum is from Laine & 
Gottesman (1998). (b) Azimuthally averaged radial profiles of the CO (275 spatial resolution) and H I (15" spatial 
resolution) emission in NGC 7479, together with the total (molecular + neutral hydrogen) gas profile. The vertical 
scale gives the mass surface density in Mq pc“^. The molecular profile is shown with filled triangles, connected 
with dashes, the H I profile with stars, connected with dots, and the total gas (molecular + neutral hydrogen) 
profile with open circles, connected with a solid line. 
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be expected if it followed the rotation curve derived 
from the gravitational potential, obtained by convert¬ 
ing the itl-band surface brightness into a mass sur¬ 
face density ( Quillen et al. 199^ . Since Quillen et al. 
(1995) did not consider the three-dimensional nature 
of the bulge, and the M/L ratio as determined from a 
iiT-band image may be variable near the nucleus, the 
gravitational potential and the rotation curve in the 
central are uncertain. 


The remarkable departures of the observed emis¬ 
sion from the expected location traced by the line in 
Fig. I outside the nuclear region indicate that the gas 
motion along the bar is inconsistent with pure circular 
rotation. Assuming that the northwestern side of the 
galaxy is the near side (spiral arms trailing) and that 
there are no substantial vertical motions, it is possible 
to make a rough estimate of the radial velocity of the 
gas along the bar, taking the tangential speed of the 
gas in the bar from the bar rotation rate. There will 
an additional component of tangential speed, but if 
the gas is captured in orbits that are elongated along 
the bar, a first order estimate of the tangential ve¬ 
locity is provided by the bar pattern speed. Using 
the bar pattern speed of 27 km s“^ kpc“^ (Laine 
et il. 199^ ), the tangential velocity at the distance 
of 15" (3.2 kpc) is 65 km s“^. The observed veloc¬ 
ity at this point in the clump north of the center is 
only slightly higher, 2370 km s“^, than the velocity 
of the nucleus. Accounting for the difference between 
the directions of the gas bar major axis and the line 
of nodes (about 31°) and the inclination (51°; Laine 
& (lottesman 1998), the observed emission velocities 
can be reproduced with a radial inflow component of 
110 km s“^. 


The plot along the CO major axis (Figure ||a) also 
gives a hint of the existence of two separate kine¬ 
matic components. The first component shows up as 
a large linear velocity gradient across the central few 
arc seconds (190 km s“^ arcsec“^ or 1.2 km s“^ pc“^) 
that is consistent with solid body rotation. After cor¬ 
recting the observed velocities for the inclination, the 
maximum rotation velocity in the nuclear disk is 245 
km s“^ at a radius of 1". The noncircular kinematic 
component can be seen as extensions toward lower 
declinations and lower velocities (marked with El in 
Fig. I a) and higher declinations and higher velocities 
(marked with E2 in Fig. ||a). 

We have also made position-velocity maps perpen¬ 
dicular to the bar (at PA 90°) at several positions 
along the bar, using both the CO and the Ha data 
cubes. The resulting maps, together with an image 
of the integrated CO and Ha emissions, are shown in 
Figure |[ The fluxes between the horizontal bound¬ 
aries indicated in the map were summed together to 
improve the signal-to-noise ratios within every num¬ 
bered region. The velocity resolutions are 10.4 km s“^ 


(CO) and 12.1 km s ^ (Ha). 

The overlay of the total CO intensity (at 2'.'5 spa¬ 
tial resolution) contours on an Ha image (at O'.'S spa¬ 
tial resolution, kindly given to us by Dr. J. Knapen) 
reveals that the CO and Ha emission peaks occur 
mostly at the same locations within the available spa¬ 
tial resolution. However, in regions 3 and 7, the Ha 
emission maxima are at a larger radius and larger 
distance from the bar major axis than the CO max¬ 
ima. The position-velocity maps show that in Ha, 
two velocity “plateaus” of emission, both upstream 
and downstream from the dust lane, are connected 
by a steep CO (and Ha) velocity gradient. Strong 
CO emission only exists in the region of the steep 
velocity gradient. 

The largest velocity gradient in the CO position- 
velocity maps presented in Fig. ||is about 240 km s“^ 
arcsec”^ in region 7. Correcting this for the inclina¬ 
tion gives a gradient of 1.9 km s“^ pc“^. The least 
steep gradients further out along the bar have values 
around 0.5 km s“^ pc“^. 

4. Discussion 

4.1. Gas flow and molecular gas in the bars 
4-1.1. Comparison of observations and models 

Observations of NGC 7479 show dust lanes dis¬ 
placed toward the leading edges of the bar, accom¬ 
panied by CO emission (Fig. which traces dense 
compressed molecular material. The large velocity 
gradients that we see in our observations occur along 
these gas/dust lanes. This is consistent with the gen¬ 
eral expectation of shocks in the compressed interstel¬ 
lar medium. In shocks the component of velocity per¬ 
pendicular to the shock has an abrupt change in mag¬ 
nitude. The large velocity gradients in our observa¬ 
tions are also consistent with the models which show 
that gas streaming along the bar has an apogalacticon 
near the bar dust lanes, where the gas flow changes 
from having a radially outward directed velocity com¬ 
ponent to having a radially inward directed compo¬ 
nent (see, e.g.. Fig. 5 in Roberts, Huntley, & van Al- 
bada 1979). The line of sight at which the bar in 
NGC 7479 is observed is intermediate between per¬ 
pendicular and parallel to the bar, and thus favorable 
for seeing the large gradients in the line of sight ve¬ 
locity across the bar dust lanes. 

Figure shows position-velocity maps from obser¬ 
vations and from a Smoothed Particle Hydrodynam¬ 
ics (SPH) bar flow model ( [Laine et al. 199"^ ). The 
position-velocity slices in both the observations and 
the models were taken across the bar at the same dis¬ 
tances from the nucleus, and are thus directly compa¬ 
rable. The model data were taken from the run with 
the best fitting bar pattern speed (27 km s“^ kpc“^). 
The two-dimensional simulation frame was projected 














Velocity (km/s) 


Fig. 6.— The major axis CO position-velocity plot. The contours are at 2, 3, 4, 6, 8, 10, 14, 18, 22, 24 times the 
rms noise in the channel maps (see Table The extensions of the emission towards the south and north have 
been marked with El and E2, respectively. A possible origin of these features is discussed in the text. 
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Fig. 7.— (a). The zeroth moment plot (Fig. ||). (6). The position-velocity plot along the bar at PA 177° drawn 
at the same scale in x-axis as (a) so that the two can be compared easily. The rotation curve derived from the 
gravitational potential obtained with the help of a iF-band image (Quillen et al. 1995) has been drawn with a 
solid line in the position-velocity plot (after correcting it for the PA and the inclination of the plot) to show the 
departures from pure circular motion. 
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Fig. 8.— An overlay of the CO emission intensity on a gray scale image of the Ha emission is shown on the left 
(the O'.'S resolution Ha image was kindly given to us by Dr Johan Knapen). The contours of the CO emission 
are given at 40, 120, 200, 280, 360, 440, 600, 760, 920, 1080, 1240, 1400, 1560, 1720 H 2 atoms cm“^, using the 
standard conversion factor between CO intensity and H 2 surface density. The regions in which the emissions were 
integrated for the position-velocity maps are separated by solid lines, and numbered from 1 to 8. The corresponding 
position-velocity maps are shown on the right, the first column giving the Ha maps, the second column the CO 
maps, and the third column both overlaid, using gray scales for the Ha emission. The velocities are heliocentric. 
The contours on the position-velocity maps are fractions of the peak intensity in a map, chosen to give a clear 
presentation of the emission. The strongest CO and Ha emission occurs in a narrow ridge with a large velocity 
gradient. Weaker Ha emission is also observed both upstream and downstream from this ridge. 
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Fig. 9.— (a). Gas position-velocity maps of region 3 in Fig. |g. The Ha emission is shown with gray scales and the 
CO emission with contours. (&). A position-velocity map of a region corresponding to (a) in the most successful 
model in Laine et al. (1998), at t = 9 x 10® yr. (c). As (a), but showing region 7 in Fig. I (d). As (6), but 
corresponding to region 7 of Fig. |[ A comparison of observations and models shows that the model has succeeded 
in reproducing the observed morphology. Since the shock in the model takes place along the dust lane, we conclude 
that the observed large velocity gradients are consistent with a strong bar shock. 
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to the orientation of NGC 7479 (inclination 51°, PA 
22°; Laine & Gottesman 199S; the angle between the 
stellar bar and the kinematical major axis in the plane 
of the galaxy 25°). The spatial resolution of the model 
was smoothed by a Gaussian to 2'.'5 (400 pc) and the 
velocity sampled at 10.4 km s“^ bins, comparable to 
the GO (and roughly, the Ha) observations. 

SPH represents total gas, and therefore it traces 
both the compressed regions of dust lanes and dense 
molecular gas as well as the more diffuse ionized gas. 
In both the models and the observations (Figures |^, 
H 1^, the observed velocity changes by about 
150-200 km s“^ across the highest gas concentra¬ 
tions (which in optical and near-infrared observations 
are seen as a dust lane; Fig. |^), but the Ha emis¬ 
sion “plateaus” upstream and downstream from the 
du st lane have relatively small velocity gradients. Our 


kin ematical resiilts support the expectation that the 

bar dust lane traces shocks in the gas flow. 

The models also show the velocity “plateaus” of gas 
upstream and downstream from the location of the 
large velocity gradient. The downstream “plateaus” 
suggest that some fraction of the gas has moved out 
of the shock region. However, the mean velocity of 
the shock regions is consistent with radial inflow mo¬ 
tion as mentioned in § 3.2. Therefore, apart from gas 
that appears to have moved downstream, away from 
the shock front, our observations are consistent with 
the picture given by Regan, Vogel, & Teuben (1997, 
RVT; see also Regan 1996) who claim that all gas that 
encounters the dust lane subsequently flows down the 
dust lane into the nuclear region. 

Our observations and the comparison to the SPH 
model have shown that the response of gas (including 
molecular) to the forces along the bar in NGG 7479 is 
more consistent with the behavior of dissipative dif¬ 
fuse gas under hydrodynamic forces (pressure, viscos¬ 
ity) than the response of giant dense molecular gas 
clouds which might be expected to react like a ballistic 
particle. Both the grid-based ideal gas hydrodynamic 
codes and SPH codes have modeled gas flows, includ¬ 
ing the compression regions or shocks, satisfactorily 
in most cases. “Sticky particle” (SP) codes smear the 
shocks over relatively large spatial scales, depending 
on the particle size, and may have difficulties mod¬ 
elling the 100 pc scale, steep velocity gradients in the 
bar dust lane of NGC 7479. It is unclear at present 
which code is the most realistic in modeling the be¬ 
havior of the interstellar medium (ISM). The equation 
of state of the ISM is unknown and the multiphase 
numerical modeling of the gas has not been done yet. 

4-1.2. Properties of molecular gas in the bars 

The star formation and gas flow properties along 
bars would be easier to understand if the physical 
state of the CO emitting regions was known. Specif¬ 


ically, we would like to know if the CO emitting 
regions are gravitationally bound. RVT claim that 
the excellent agreement between their observations of 
NGC 1530 and the hydrodynamic models based on 
diffuse ideal gas suggests that the gas they traced with 
Ha kinematics is diffuse gas and not concentrated in 
gravitationally bound giant molecular clouds (GMCs; 
masses up to a few times 10® Mq) or even larger and 
more massive GMAs, which might behave more like 
ballistic particles than diffuse gas. Our observations 
and the comparison to models also suggests that the 
molecular gas in the bar of NGC 7479 does not behave 
like ballistic particles. In the disk of spiral galax¬ 
ies the situation appears different, as complexes of 
GMAs have been found in at least the nearby galaxies 
M 51 (Rand 1993a, 1993b; Vogel, Kulkarni, & Scoville 
1988) , M 100 (iRand 19951) , and NGC 4414 (|Sakamoto 
19961) . 


The CO luminosities (masses) and diameters of the 
CO emission regions along the bar in NGC 7479 are 
roughly similar to those seen in the disks of other 
spirals such as M 51, but the kinematics along the bar 
of NGC 7479 are very different. NGC 7479 has very 
large velocity gradients of 1.9 km s“^ pc“^ across the 
emitting regions. Therefore, the CO emission regions 
along the bar are unlikely to be self-gravitating. 

The virial masses of the CO emission regions in 
Mq can be roughly estimated from 


Mvir = 550i:>(aid)^ 


( 2 ) 


( Rand 1993a ; Scoville fc Sanders 1987 ). Here D is 
the observed FWHM diameter in parsecs, and aid 
is the one-dimensional velocity dispersion in km s“^. 
This formula assumes a 1/r density profile for the 
gas within a GMA. The mass required for the GMA 
to be self-gravitating is half of this mass. Using a 
velocity dispersion of 150 km s“^ and a FWHM di¬ 
ameter of 150 pc (about 1"), the resulting virial mass 
is about 2 x 10® Mq. For comparison, the observed 
masses of the clumps along the bar (assuming the 
standard conversion factor) range from 9 x 10® Mq 
to 2x 10® Mq. Thus it is unlikely that the CO gas 
concentrations that lie along the bar are virialized or 
even gravitationally bound. 


4.2. Star formation in the bar of NGC 7479 

Star formation rates (SFRs) from the Ha fluxes in 
the offset gas and dust lane along the bar have been 
calculated from 


SFRna = 7.07 x 10"'‘^L(Ha) Mq yr 




(3) 


( Bushouse 1987 see also Kennicutt 1983 ). Equation 
(3) assumes a Salpeter initial mass function (IMF) 
down to stars of mass 0.1 Mq and up to stars with 
masses of 100 Mq. The Ha fluxes have not been 
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corrected for extinction, so the SFRs from equation 
(3) are lower limits. The integrated SFR in the bar 
region is around 0.2 Mq yr“^. Uncertainties due to 
extinction and the choice of the IMF with its cut¬ 
offs are discussed in Bushouse (1987). Specifically, 
Bushouse compared the global star formation rates 
as calculated from the Ha fluxes (uncorrected for ex¬ 
tinction effects) to global far-infrared fluxes. The lat¬ 
ter were found to be three times larger on average, 
and therefore an average extinction correction factor 
of three applies to isolated galaxies. For comparison, 
Martin & Friedli (1997) measured the ratio of the Ha 
and H/3 lines in a few selected H II regions along the 
bar, and obtained a bar extinction value of 0.8 mag in 
Ha, which corresponds to a factor of 2. Multiplying 
the measured Ha flux by 2.5 results in a star forma¬ 
tion rate [SFR(Ha)] of 0.5 Mq yr“^ along the bar of 
NGC 7479, comparable to the value (0.4 Mq yr“^) 
obtained by Martin & Friedli (1997). The SFR for the 
whole galaxy, calculated from the infrared luminosity 
[L(IR); 8.9x 10^° L(7); [Young et al. 1989|, is about 
45 Mq yr-|^ [SFR(IR) = 1.34x10"^^ T(IR); [Hunter 
et 4l- 1981 ; Bushouse 1987], 90 times larger than the 
SFR along the bar and thus largely not seen in op¬ 
tical tracers, although the Ha luminosity of the bar 
region is only about 10% of the total Ha luminosity 
of NGC 7479. 

Bushouse (1987) finds that SFR(IR) is 3-6 times 
larger than SFR(Ha) among normal galaxies, (3 times 
in isolated, 6 times in interacting galaxies), most 
likely because of extinction in the optical. Yun & Hib¬ 
bard (1998) confirmed the significantly larger SFR(IR) 
values and also found that the Ha component seen 
in the nuclear starburst systems is dominated by the 
starburst superwind and scattered emission. It is pos¬ 
sible that the nucleus of NGC 7479 hosts a starburst 
behind dozens of magnitudes of optical extinction. 
About 15% of the 21 cm continuum emission is con¬ 
centrated in the central beam (4'.'3 x 3785), giving 


a r ough size scal e for the starburst region (jLaine fc 
Go desman 1998 ). 


The large velocity gradient of 0.5-1.9 km s“^ pc“^ 
across the bar dust lane of NGC 7479 is likely to 
be dynamically important and affect star formation. 
The large velocity gradient can prevent large gas con¬ 
centrations within the bar dust lane from being self- 
gravitating, as indicated by our simple calculation in 
§ 4.1.2. 

The kinematical conditions in the bar dust lane are 
quite different from those in the spiral arms of M 51, 
which has unusually strong spiral density waves with 
small pitch angles. Although there is a density wave 
induced velocity gradient of 0.05-0.1 km s“^ pc 


^-1 


across the front half of the spiral arms of M 51 (Vogel 
et al. 1988; Rand 1993a), this gradient is an order of 
magnitude smaller than that in the bar of NGC 7479, 


and is furthermore in a direction which reduces the 
amount of shear which otherwise occurs in a disk with 
a flat rotation curve. 

Thus, in spiral arms the largest concentration of 
dense gas and the regions of reduced shear occur in 
the same places, and both effects increase the suscep¬ 
tibility of gas to gravitational collapse. In bars the 
largest concentration of dense gas occurs in a region 
of a large velocity gradient, and this opposes the sus¬ 
ceptibility of gas to gravitational collapse. 

This could be the reason why star formation is ob¬ 
served upstream and downstream from the dust lane 
as well as within it. Although the regions upstream 
and downstream from the shockfront have lower gas 
density, their smaller velocity gradients may enable 
star formation. As discussed in § 3, the regions up¬ 
stream and downstream may have molecular gas that 
has a column density which is too low to be detected 
when averaged over the synthesized beam, but the 
density may be high enough to enable star formation, 
making it easier to understand the wide distribution 
of Ha emission across the bar. 

The regions with the highest gas density are re¬ 
gions with the largest velocity gradients. Although 
stars are certainly forming in this dense dust lane gas, 
the star formation rate may be limited by the large 
velocity gradients: the long timescale for star forma¬ 
tion along the bar (2 x 10® Mq/0.5 Mq yr“^ = 
4.0 X 10® yr) implies that the bar gas is not experi¬ 
encing a starburst. Again, this conclusion is subject 
to the effects of extinction, as discussed before, al¬ 
though these effects are unlikely to change the basic 
conclusion. Elmegreen (1979) and Kenney & Lord 
(1991) have also suggested that the tidal forces on 
the gas clouds along the bar dust lanes may inhibit 
star formation. 


4.3. Nuclear gas mass fraction 


The molecular gas mass in the centers of some spi¬ 
ral galaxies is a large fraction (up to at least 50%) 
of the total dynamical mass ( [Turner 1994 ; Kenneyj 
1997). Under these conditions, the self-gravity of 


the molecular gas becomes important and may drive 
further inflow towards the center (^hlosman et al 


1118S;_ShlosmaiQg1^1. ljl90yElgl]er_fc_^hlo sman 1994 


Knapen et al. 1995 ; Wada fc Habe 1992| ). The gas 
mass fraction in the central disk of NGC 7479 can 
be calculated in the region that exhibits circular ro¬ 
tation, using the standard CO luminosity to H 2 mass 
conversion factor. The rotation velocity in the cen¬ 
ter is about 245 km s“^. We use a radius of 1" for 
the nuclear gas disk. Using a distance of 32 Mpc, the 
total dynamical central mass becomes 


Mdyn =2.2 X 10®Mq. 

Gr 


( 4 ) 
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The molecular mass in this area is 3.3 x 10® Mq 
or 15% of the dynamical mass. The molecular gas 
mass is a considerable part of the total mass in the 
center, although it does not dominate the dynamics. 
If the rest of the mass is stellar, the average mass- 
to-luminosity ratio in the if-band in this central area 
is 0.7 in solar units. This would imply a relatively 
young (a few Gyrs) stellar po pulation or alte rnatively 
a metal rich old population (Worthey 1994). 


Since the star formation rate along the bar is 
0.5 Mq yr~ ^ and the gas infl o w rate is a t least 
4 Mq yr“^ ([Quillen et al. 19^ ; Laine 1996 ), most 
of the gas presently in the bar region should reach 
the circumnuclear region before it turns into stars. 
The uncertainties are caused by the unknown amount 
of extinction, which may severely suppress the ob¬ 
served Ha fluxes along the bar dust lane and cause 
an underestimate of the bar SFR. The net inflow of 
2 X 10® Mq of gas would double the nuclear mass 
in less than a Gyr, meaning that NGG 7479 is in a 
short-lived evolutionary state. The resulting increase 
in the nuclear mass would likely drive the evolution 
of the galaxy towards an earlier Hubble type (Kor- 
mejidy 1993|; [Friedli & Benz 1993|; [Kenney & Jogee 
19S|^^ 


5. Conclusions 


The high spatial resolution GO observations of 
NGG 7479 together with Ha data have shown the 
gas distribution and motions in a dynamically young 
barred spiral g alaxy (cf. [Kenney fc Jogee 1997 ; Laine 
Heller 1998 ). Our main conclusions are summa¬ 
rized in the following: 


1. The GO emission has a strong central peak, but 
about 50% of the emission comes from regions outside 
the nuclear zone, mostly from a clumpy distribution 
along the bar dust lane. 

2. The bar dust lane and the GO gas emission 
occur at the same locations within the resolution of 
the observations. 


3. The GO velocity field is consistent with circular 
rotation only in the central 2" (300 pc) diameter re¬ 
gion. Outside this region, the velocity field deviates 
strongly from pure rotation and the velocities must 
have a substantial (at least 100 km s“^) radial com¬ 
ponent along the bar. 

4. The GO emission across the bar dust lane has a 
large velocity gradient (up to 1.9 km s“^ pc“^) with a 
projected velocity change up to 200 km s“^. The Ha 
emission shows a similar large gradient, but it also 
has “plateaus” of constant velocity emission both up¬ 
stream and downstream from the large velocity gra¬ 
dient. The large velocity gradient is consistent with 
combined effects of a shock and a change in the radial 
direction of gas streamlines in the bar dust lane. 


5. The gas clumps along the bar have masses from 
9x10® Mq to 2x10® Mq if the standard GO-to- 
H 2 conversion factor is assumed. However, the large 
velocity gradient across these complexes in the bar 
dust lane makes this assumption suspect. 

6. No strong GO emission exists in the bar region 
outside the dust lane, and the gas complexes in the 
dust lane are unlikely to be gravitationally bound be¬ 
cause of the large velocity gradients across them. 

7. Based on the relative strengths of the GO and 
Ha emissions in the nucleus and along the bar, we 
conclude that either the optical emission from the 
nucleus is strongly reduced by extinction or the star 
formation efficiency in the nucleus is lower than in 
regions along the bar. 

8. The star formation timescale along the bar is 
long, a few Gyr. The large velocity gradients in the 
bar dust lane may prevent the bar gas from forming 
stars more rapidly. 

9. Along the bar, the gas inflow rate of 4 Mq yr“^ 
is much greater than the star formation rate of 0.5 Mq 
yr“^, implying that most of the gas presently in the 
bar will reach the circumnuclear region. If that hap¬ 
pens, within the next ^ 1 Gyr the central mass of 
NGG 7479 should roughly double. However, extinc¬ 
tion in the bar dust lanes and in the nuclear region 
causes the cited SF rate to be a lower limit, and the 
actual amount of gas reaching the circumnuclear re¬ 
gion is likely to lower. 

10. The nuclear gas mass fraction is likely to be at 
least 15% of the total dynamical mass within the cen¬ 
tral gas disk, and therefore, the molecular gas affects 
the dynamics in the nuclear disk. 
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Fig. 10.— The contours of the naturally weighted 
CO emission overlaid on the gray scale image of a 
J — K color map (Laine 1996), using data taken at the 
United Kingdom Infrared Telescope (UKIRT). Darker 
colors are redder regions in the gray scale map. The 
contour levels are at (4, 8, 16, 32, 48, 96, 144) x 10^^ 
H 2 molecules cm“^. The excellent spatial correlation 
shows that strong CO emission is associated with the 
high extinction in the bar dust lane. 
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